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Abstract 

It has recently been pointed out that under certain circumstances the back-reaction of 
charged, massive Dirac fermions causes important modifications to AdS2 spacetimes arising 
as the near horizon geometry of extremal black holes. In a WKB approximation, the modified 
geometry becomes a non-relativistic Lifshitz spacetime. In three dimensions, it is known that 
integrating out charged, massive fermions gives rise to gravitational and Maxwell Chern-Simons 
terms. We show that Schrodinger (warped AdSs) spacetimes exist as solutions to a gravitational 
and Maxwell Chern-Simons theory with a cosmological constant. Motivated by this, we look for 
warped AdSs or Schrodinger metrics as exact solutions to a fully back-reacted theory containing 
Dirac fermions in three and four dimensions. We work out the dynamical exponent in terms of 
the fermion mass and generalize this result to arbitrary dimensions. 



*jhung@perinieterinstitute.ca 

^djatkar@perimeterinstitute.ca 

■f-asinhaQperimeterinstitute.ca 



1 



Contents 



1 Introduction 2 

2 Schrodinger spacetimes in Chern-Simons theories 4 

3 Schrodinger spacetimes from fermions 5 

3.1 2+1 Dimensions 7 

3.2 3+1 Dimensions and beyond 11 

3.2.1 Massless fermions 13 

3.2.2 Massive fermions 13 

3.3 Can there be an exact Lifshitz solution? 14 

4 Discussion 15 

A Conventions 17 

B Spin connection and the Dirac Equation 18 

C Stress tensor of fermionic probes in AdS 19 

C.l AdS2 xR 19 

C.2 AdS2 X i?2 20 

C.3 AdS2 xS^ 21 

D A toy model justification for our semi-classical treatment 21 



1 Introduction 

In recent times, there has been a concerted effort to find real world apphcations of the AdS/CFT 
correspondence. By extending the rules of the correspondence to finite temperature, there are 
now well established prescriptions to calculate retarded Green functions which enables us to 
extract transport properties of interesting physical systems. Although none of these systems, 
strictly speaking, describe the real world, it is hoped that by studying these problems, one 
will gain intuition about real physical systems which are for example in the same universality 
class. On one hand, AdS/CFT methods have found apphcations to particle physics, e.g., the 
viscosity of the strongly coupled quark gluon plasma [1]. On the other, these methods have been 
extended to studying table-top condensed matter physics [2]. It is fair to say that applications 
of gauge/gravity duality will yield many more surprises. 

One application that has garnered much attention lately is the computation of fermion Green 
functions using AdS/CFT which yield non-Fermi liquid type behaviour [3]. In the simplest 
example, one starts with a Reissner-Nordstrom AdS black hole and considers the Dirac equation 
in this background. In the zero temperature limit, the black hole becomes extremal and the 
near horizon geometry contains an AdS2 factor. This AdS2 is thought to play an important 
role in the existence of a Fermi surface. The resulting fermion Green function shows non-Fermi 
liquid type behaviour which is exciting in the context of connections with condensed matter 
systems, e.g. high Tc superconductivity. However, in the zero temperature limit, the black hole 
and hence the corresponding dual field theory, still has a large entropy which is far from what 
is expected in the realistic situations. 

Recently in [4], it has been pointed out, that in certain regimes the back reaction of the 
fermions cannot be ignored. In fact, the AdS2 geometry has been argued to be modified into a 
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"Lifshitz" [5] type and the associated entropy enigma is thought to disappear since the horizon 
area is now zero. Motivated by this, we examine the problem of back-reaction of Dirac fermions 
in 2+1 and 3+1 dimensional gravitational theories in the context of AdS/CFT in some detail. 
The 2+1 dimensional bulk theory will describe 1+1 dimensional field theories. In this context 
interesting non-Luttinger type behaviour has been found in [6]. The 3+1 dimensional bulk the- 
ory will describe interesting 2+1 dimensional field theories which are of relevance to condensed 
matter physics related to high Tc superconductivity. Wc will restrict our attention to the sim- 
plest case where only the metric, a U (1) gauge field and a charged, massive Dirac fermion is 
involved. 

We begin with the discussion of a 2+1 dimensional bulk theory. It is known that [7, 8, 9] 
when one integrates out a charged fermion in 2+1 dimensions, the effective action that one 
is led to contains the Maxwell and gravitational Chern-Simons terms. Gravitational Chern- 
Simons theories have been the subject of much attention recently due to the connection with 
topologically massive gravity (TMG)[10]. In TMG null warped AdSa have already been found. 
The new finding in our case is their existence in the presence of U{1) Chcrn-Simons term 
A A F. We show that in addition to the usual asymptotically AdSs spacetime, one also gets 
Schrodinger spacetimes [11] whose dynamical exponent depends on the Maxwell Chern-Simons 
or the gravitational Chern-Simons terms. Similar solutions (called warped AdSa) have already 
been reported in [12] in the presence of the gravitational Chern-Simons terms. Inspired by this 
finding, we look for Schrodinger spacetimes as exact solutions to back-reacted Dirac fermions. 
Since fermions obey the Pauli exclusion principle, they cannot be treated classically. We will 
treat the system of gravity, gauge field and fermions in a semi-classical manner following [13, 14]. 
In this approach, the fermion stress tensor and current appearing in the equations of motion are 
evaluated as expectation values in some state We will consider \Q) to be the state made 

from N fermions. It will turn out that the role of h is played by 1/Nq so that for a fixed charge, 
taking the large limit will correspond to making h small which is needed for a semi-classical 
approximation to make sense. The solution to the full problem in a self-consistent manner in 
general is a very hard task. Typically one resorts to some approximation as in [14] where a 
WKB approximation is used to derive the Oppenheimer-VolkofF equations. Quite remarkably, 
our problem will turn out to be amenable to an exact solution! The key ingredient that makes 
this possible is the fact that the fermion couples to the gauge field. The equations of motion 
then lead to certain important constraints which make an exact solution possible. In arbitrary 
bulk dimensions, D, we will show that the dynamical exponent is given by 

z = ±2mL - (L> - 2) , (1.1) 

where L is the length scale entering in the definition of the cosmological constant in the usual 
way. Lifshitz metrics are known to emerge as a WKB approximation [4] . It is not known if they 
are exact solutions to a system of self-gravitating fermions as in the case of the Schrodinger 
metrics. In the 2+1 dimensional case, as we will show, Lifshitz metrics are not exact solutions 
to the equations of motion. However, in higher dimensions wc expect there to be more general 
ways of solving the equations and it could be that other interesting solutions exist. We leave 
this interesting question as an open problem.^ 

This paper is organized as follows. In section 2, we consider Chern-Simons theories in 2+1 
dimensions and show that Schrodinger type metrics with specific dynamical exponents exist 
as solutions to the system of equations. Motivated by this, we ask if similar metrics exist 
as fully back-reacted solutions to Dirac fermions in section 3. We work out the dynamical 

^This method is used quite routinely in nuclear physics. See, e.g., [15]. 

^The fermionic stress tensor eq.(7.43) in [4] can be shown not to lead to Schrodinger spacetimes. 
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exponent in arbitrary dimensions, i.e., where the original metric would have been AdS2X R^~^. 
We conclude with open questions in section 4. In appendix A we spell out our conventions 
and in Appendix B, we derive spin connections for a certain class of metrics, which includes, 
AdS, Schrodinger and Lifshitz metrics and set up the Dirac equation in these backgrounds. In 
Appendix C we consider the stress tensor for a Dirac fermion^ in a background containing an 
AdS2 factor. We denote the charge of the fermion by q and mass by m. The AdS2 arises as 
the near horizon geometry of some extremal Reissner-Nordstrom black hole. In particular, a 
gauge field At needs to be turned on for Einstein's equations to be satisfied. We consider a 
charged, massive Dirac fermion in this background as a probe and work out the stress tensor. 
We find that quite generally an off-diagonal component of the stress-tensor is turned on. We 
show that when m < \Jl(f- + 1, the back-reaction of fermions becomes important. In appendix 
D, we study a quantum mechanical toy model in path integral formalism which corresponds 
to coupling of macroscopic fermi gas to a bosonic harmonic oscillator. We integrate out the 
fermions and obtain an effective action for the boson. We show that when there are a large 
number of fermions, the back reaction of fermions is well approximated by replacing fermion 
number operators by their expectation values evaluated in the N fermion state. 



2 Schrodinger spacetimes in Chern-Simons theories 

We will begin with the discussion of dimensional Einstein-Maxwell-Chern-Simons theory. 
Our motivation for looking at this theory came from the observation that for a probe fermion 
in AdS2 geometry, off-diagonal components of the stress tensor are non-zero. (This has been 
discussed in detail in appendix C.) As a result it may be expected that any non-relativistic 
modification of the metric will be of the Schrodinger type rather than Lifshitz. As explained 
in the introduction, on integrating out fermions in 2-|-l dimensions, one gets Maxwell and/or 
gravitational Chern-Simons theories. Rather than looking at the more complicated problem of 
back-reaction of fermions, we wish to begin by investigating if Schrodinger type solutions exist 
for an Maxwell and/or gravitational Chern-Simons theories. Consider the action given in eq. 
(A.l). The equations of motion for the metric read 

111 1 

- -j9a{iR - J29afi - -{Fa^F^^ - -F'^gaji) + IJ-cCa/i = , (2.1) 

where C^*^ = e^^'^^V a{R^p — is the Cotton tensor. The equations of motion for the gauge 

field read 

V^F^'i' + 2^ipe^''"F,^ = Q. (2.2) 
It can be shown that the following Schrodinger^ (or warped AdSa [12]) solutions exist 

ds'^ = -r''dt^TArdtdx + , (2.3) 

4r^ 

At = ^[rM^-l)(l^/XG±2/XG-^)]'/^ (2.4) 

z = 0,1,t4mf,^^|^, (2.5) 



^Dirac systems in AdS spacetimes have been studied also in [16]. 

^This nomenclature is motivated by the obvious similarity of the metric with the standard Schrodinger spacetimes 
which involve additional spatial directions. Strictly speaking this is an abuse of nomenclature since in our case there 
are no spatial directions which transform under dilatations. The conformal boundary of such spacetimes has been 
discussed in detail in [17]. 



4 



To avoid clutter, we have set L = 1. Here z is the standard dynamical exponent and z = \ 
corresponds to the AdS solution, whereas z = corresponds to a chiral wave AdS solution. 
The thing to note here is the existence of Schrodinger solutions when either or both of the 
Chern Simons terms exist. The existence of these solutions in the presence of the gravitational 
Chern-Simons term can be anticipated from the general discussion in [18]. The case z = 2, 
(^^ = ::pl/2 or /XG = =Fl/3) is called the null- Warped AdSa [12]. 

At this point it is worth noting that Chern-Simons terms are not invariant under parity (P) 
and time reversal (T) transformations but are invariant under combined transformation PT. This 
symmetry is also shared by the Schrodinger background, which under P or T transformation 
swaps the sign of dtdx term but the metric is invariant under combined PT transformation. 
However, in case of z = 1 the metric is diffeomorphic to AdS metric and preserves P and T 
separately. This does not contradict earlier conclusion because for 2 = 1, all components of 
gauge field vanish and in that case it is natural to expect to recover P and T symmetry. 

With a bit more work it can also be shown that a Lifshitz metric is not a solution to the 
above set of equations. We begin with the ansatz 

d^' = -^ + At = (P{r),Ar = 0,A, = x{r). (2.6) 

Then the tt component of the metric equations of motion leads to r^"*"^^!^'^ + r^'x'^ = while 
the XX component leads to 4 — -|- r^+^^c/)'^ -|- r*^x'^ = 0. Combining these two we are led 
to z = ±1. The solution z = 1 \s the usual AdS. The choice z = —1 can be shown to lead to 
imaginary gauge fields. Thus we conclude that Lifshitz is not a solution to this system. We now 
turn to the more complicated problem of considering back-reaction due to fermions. 



3 Schrodinger spacetimes from fermions 

If we compute the stress tensor for a probe fermion in an AdS2 background as in appendix C, 
we find that it contains off-diagonal components. This suggests that if their back-reaction is 
taken into account an exact background geometry of the form of a Lifshitz geometry, which is 
diagonal, is unlikely. Given however, the intuition from flat space that fermions and the Chern- 
Simons terms are intimately related, and the similarity in structure of the energy-momentum 
tensors of the gauge fields in the presence of the gauge Chern-Simons term and that of the 
fermions, one is tempted to conjecture that the back-reaction of charged fermions would also 
lead to Schrodinger spacetimes, exactly as shown to happen, in the previous section when the 
back-reaction of gauge fields in the presence of Chern-Simons terms is taken into account. To 
be explicit, the full Lagrangian we will be considering is (total bulk spacetime dimensions is 
denoted hy D = d+1) 



S = [ d'^+^xJ^ 



R-2A- ^Ff.^Fi"' + iiipVei^D^iP - mipip) 



(3.1) 



where, A = — '^^2L'^^ ' P-matrix convention is given in Appendix B. Note here that in general 
the fermion field will be a sum over modes 

ip{r,t,x) = ^ {akipa,k{r,t,x) + blil;b,k{r,t,x)j , il;\r,t,x) = ^ (4V'l,fe ('''*' a;) -|- 6feV'6,^(r, t, x)) 

k k ' ' 

(3.2) 

where ak,a\,hk,b^i^ denotes the anti-commuting creation and annihilation operators of the cor- 
responding fermion and anti-fermion modes ipa/b,k{''^jt,x) respectively, k denoting some general 
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quantum numbers of the modes, each mode i^a/b,k {^i ^i ^) ^ two-component spinor being a so- 
lution to the Dirac equation, and that to lowest order ignoring quantum fluctuations of the 
background we have, 

' 4, } = ^^h,k, = {h, , 4, }' (3-3) 
while all other anti-commutations between operators vanish. Since the background geometry 
itself is strongly back-reacted by the fermions, in general it is slightly ambiguous to define what 
the vacuum would be. We assume the existence of such a stable vacuum such that 

ak\0)=bk\0) = 0. (3.4) 

Any general excited state would be of the form 

|Q),enera;=n4.n^l|0)' M 
i j 

although they are not generally the AT-fermion ground state at zero temperature, and are thus 
not necessarily stable against decay. The Einstein equation is, in our conventions (Appendix 
A), 

Ra^ - ^gaCsR - ^^^^9aP ' ^(^a7 V " i^'^a^) = ^(^I^^IQ) ' (3-6) 

where is as defined in (C.3) and we reproduce here for convenience 

= \ {-ii>T^^D,^i^ + #^(^r,)^) . (3.7) 

This form of the stress tensor is given in [19] and to the best of our knowledge was first given 
in [13]. The covariant derivative Dfj_ is defined in equation (C.2). The Maxwell's equations are 
given by 

V^F^. = {Q\j^Q) = -gV^r^^. (3.8) 

Note that in the above expressions, the fermion stress tensor and electric currents that 
back-react on the background are really the expectation values of the corresponding operators, 
evaluated on the state \Q). At zero temperature, the ground state of N fermions would cor- 
respond to piling up the various modes, starting from the bottom mode. Naively, both the 
currents and the fermion stress tensor would take roughly the form J2k C'k'^k\'4'k{i^j ^) ^)P > whose 
expectation value in the state \Q), would involve a sum over contributions of modes that are 
occupied in \Q). However, here we should note that in general one is confronted with a diverg- 
ing expectation value due to contributions even of unfilled modes. It is ambiguous to define 
normal-ordered operators in curved space, since vacuum energy gravitates as well [19]. Here we 
approach the problem via analytic continuation [19]. As wc will see in the following section, the 
modes that are consistent with the background symmetry and the Gauss's law constraints are 
simply chiral plane waves along t in Schrodinger space, with trivial radial dependence. There- 
fore they contribute in the expectation value of the current or stress tensor via {Q\ J2i(4^i\Q) 
simply as 2Zf'(l) = N and thus precisely in the same way the unoccupied anti-fermion modes, 
which we have infinitely many of them will contribute via 

{Q\ b4\Q) = 1 + 1 + 1... = C(0) = (3.9) 

i 

where in the second equality we attempt to regularize it via analytic continuation. Adopting 
this approach would correspond to shifting N by — ^. Note that this is sub- leading in N, and 
in the large N limit, which is the main focus of this paper, we will ignore this term. 

The analysis that we present below is very similar in spirit to that in [14] for a system of 
self-gravitating fermions. The essential steps in this are as follows: 
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1. We want to find a particular set of fields involving the metric, gauge field and fermions 
such that it is an exact solution to the equations of motion given above. 

2. In [14], a certain ansatz for the metric compatible with some pre-supposed symmetries is 
taken. The Dirac equation is solved in this background for a complete set of modes. 

3. With these solutions, one works out the stress tensor entering Einstein's equations, 
taking care of the antisymmetric nature of the fermions. 

4. The full set of equations is now solved self consistently. 

In the regime of large quantum numbers and using a WKB approximation [14] recovers the 
Oppenheimer-Volkoff equations^. In [4], a WKB approximation was used to obtain a Lifshitz 
solution. We will not make any such approximations in what we do. In addition to the metric 
and fermions as in [14] we also have a gauge field which makes the calculation somewhat more 
complicated. To keep life simple, we will begin with the ansatz for a Schrodinger spacetime 
with an undetermined dynamical exponent. Then we will ask if there exist self-consistent set of 
solutions. To be compatible with the symmetries of the problem, we will assume that the gauge 
field is only dependent on the radial coordinate. At this stage, we note that the stress tensor 
for fermions has an h [14]. We will find that 1/Nq plays the role of h in what we do. A quick 
way to see that this is true is that in our solution, the stress tensor from fermions is of the same 
order as that of the Maxwell fields and the Einstein tensor, not surprisingly since we found a 
fully back-reacted solution. However, the general form of the charged fermion stress tensor is 
roug hly T-f ~ hNq, using the anti-commutation relation of the operators. Therefore our solution 
naturally forces Nq ~ 1/^. It is interesting to note here that for small Nq, it would imply that 
h is large and that a classical approximation would not stand. This conforms with our intuition 
that a classical approximation improves when the fermi gas becomes macroscopic. To find 
further support for our approach, we present a path-integral calculation of a toy-model consisting 
of fermions coupled to bosons in Appendix D. We found that quite generically in the limit of 
large fermion number N, our semi-classical treatment adopted here is a good approximation. 

3.1 2+1 Dimensions 

Like the Chern- Simons terms, massive fermions in 2-1-1 dimensions also break P and T symmetry. 

We therefore expect the result here to be similar to what we got in the presence of Chern-Simons 
terms. Anticipating that the back-reacted exact geometry is going to be given by the Schrodinger 
metric, we start with the ansatz metric exactly as in the previous section, namely 



where e = ±1, and leave the gauge fields and the index z to be determined by the equations of 
motion. The three sets of equations involved are the Maxwell equations, the Einstein equations 
and the Dirac equations, to be solved simultaneously. We will, unless otherwise specified, set 
L = 1 to avoid clutter in the equations. 

To begin with, we compute the spin-connection necessary to build the Dirac equations. 
Choosing explicitly the vierbein as 



ds^ = L^{-r^dt^ - Aerdtdx + —^), 



(3.10) 



e° = r~^{r''dt + 2erdx), 



= 2er^ 2 dx, 




dr, 



(3.11) 



^In [20], a similar analysis was used to study holographic neutron stars. 
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the corresponding spin-connection is 

cjOi = i^dr, oj^'^ = r-^{r''zdt + 2erdx), u^"^ = r'^ir'^iz - l)dt + 2erdx). (3.12) 
Considering a particular mode of the Dirac spinor, 

^.,.(r,t,x) = e-"(^^+'-'j;| (3.13) 

the Dirac equation for a fermion of mass m takes the form 

-(2 + AiqrAr{r))i;_^^^k{r) + 2ir-'/\u! - ^^t(r)){V'-,a;,fc(r) - V'+,a;,fc(r)} 

+ [1 + 2m - z + iA-^+'/^\k - qA,{r)M'+^^^kir) - 4rV^,^,fc(r) = 0, (3.14) 

-2ir-'/^{uj - gAt(r)){V'_,,.,fe(r) - ^+,^,fc(r)} + (2 + 4igrA(r))^+,^,fc(r) 

-[1 + 2m - z - ir^-^+'/^\k - 9^x(r))]V-,c.,fc(r) + 4rV;,^,fc(r) = 0. (3.15) 

It is simplest to start examining the Maxwell's equations (3.8). We shall begin by picking 

the gauge 

Ar = 0. (3.16) 

We will be looking at radially symmetric solution such that the metric and the gauge fields con- 
tain only r dependence. Given our gauge choice, the Maxwell's equation for the Ar component 
becomes purely a constraint on the vanishing of the r-component of the fermion current. That 
gives, explicitly, 

Jr = -Q^ ' ^^^^^ 

where we have defined the operators to be distinguished from the modes V'±,w,fe 

V' = ( 1+ ) , V'i = E a^e-^^^^'-'^^r±,.,k + ble'^^^'-'''m^,k- (3-18) 

At the same time, to preserve the symmetry of the system, we would like to keep the magnetic 
field zero i.e. F^r = 0, we allow only for 

A.,{r) = af, (3.19) 
for some constant af, which, from the gauge equations, further give 

A^(r) + rA^(r) = -qe{{r-^/'){Q\{4,l_ - ^^)(V^_ - V'+))|Q) = 0. (3.20) 
To satisfy both requirements (3.17,3.20), in general we have 

= r+,u.,ki^), V'^,.,fc(r) = ijl,^,kir), (3.21) 

for each contributing mode occupied in state vector \Q). 

Now going back to the Dirac equations, given the relation (3.21) the Dirac equations reduce 

to 

r((l - 2m + z)r+l,k + ^rV'+lfc') ± ^ '?^-)Vvi,fc = 0. (3.22) 

It is important to note that with the relation (3.21), the At and timc-dcrivativcs have 
completely dropped out from the Dirac equation. This implies that all the modes with different 
oj share the same radial wave-function, which also means 

r±,o.,k = A-,k- (3.23) 
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Also, the term proportional to {k — qA^) renders the two equations inconsistent. To obtain 
non-trivial solutions, it implies 

af = -, (3.24) 
Q 

which cannot be satisfied generally for an arbitrary A;- mode in a given background of some af. 
This implies that the modes consistent with the background are those with arbitrary u but 

fixed k such that (3.24) is satisfied. These solutions however are gauge equivalent to having 
k = Ax = 0, via the gauge transformation X = ~^ e^*'^^V) such that correspondingly 

Ax Ax — ^. We will thus from now on set Ax = k = 0. i.e., in full generality, we have 

^(^'0 = ( 1 ) EK»«"'"^**+(0 + ^i.e-'**+(0), (3.25) 

where we have 

r±,u.,o = '^±,u>,o = ^+{r)- (3.26) 

To avoid complications, for the moment we treat the spectrum as discrete. This does not 

alter our results in any significant way even when the continuous limit is taken. We would 
like to consider the back-reaction of a large number of fermions. This can be done by building 
explicitly a fermi gas by constructing the A/"-fermion state vector \Q) 

N 

IQ) = n4jO). (3.27) 

i 

Let us emphasize here that is the number of modes excited, which is infinity in a continuum 
limit for finite fcrmi-cncrgy. The fcrmion density distribution p however is smooth. As we have 
demonstrated above, since modes of different uj have the same r-wave function, evaluating the 
expectation value of the stress tensor and the current on would be proportional to the term, 
schematically given by 

N oc 1 

(Q|^iV;+|Q) = \^+(.r)\\Q\Y,ala^, + ^blbM = HN - -)\^ +{r)\' , (3.28) 

i i 

where the unfiUed-anti-fermionic states contribution bb^ has been regularized as discussed in 
(3.9). Since they are sub- leading in N after regularization, we will not include them explicitly 
in the rest of our discussion. Note that in the expectation values of the operators, cross terms 
occurring in the product i/'^Vj vanishes, and as a result both the stress tensor and the currents 
have no x, t dependence. Returning to the Dirac equation, one can then readily solve it to give 

*+(r)=piA, 5p = ^i-l + 2m- z), (3.29) 

for some suitable overall constant normalization pi^. The Maxwell equations is left with one 
component determining At- The expectation of the t-component of the current is given by 

{Q\jt\Q) = -qrHQm-f + \4>+f)\Q) = -2rt+2^-iVgn|pip. (3.30) 



^Note that because of the way we defined the measure of the frequency sum in (3.25), the normalization pi should be 
uj independent to be consistent with the anti-commutation relations (3.3) and subsequently canonical anti-commutation 
relation between the tjj and "ijj. 
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Substituting in the Maxwell equations we have 

nN\pifqA^^+^^p^ - 2rA[ - 2r^A'l = 0, (3.31) 

which gives, 

At = + 4 + c\ log r, (3.32) 

for some constants c*. The log term would appear in the energy-momentum tensor as a lone 
term. We therefore set 

c\ = 0. (3.33) 

Returning to the Einstein equations, we inspect specifically the tt component. Evaluating the 
LHS of eq. (3.6), we get 

Rtt - \guR - gtt - \{Ft^F^^ - ^F^^^) = 2(r^(-l + z)z - r^Af), (3.34) 

and the expectation value of the W-component of the fermionic energy-momentum tensor eval- 
uated on the state | Q) , 

{Q\tI\Q) = -2rH{Q\i^+dti^+\Q) - qAt{Q\i,+i;+\Q)) 

N N 

= -2ri+2^^n(^bi|V-g^tEbil') 

i i 

= -2rf+^^^bPME^^ - g^(4 - ^ + 1;,)^ ^'-''^ 

which, together with (3.32), arc substituted in to the Einstein equation (3.6). The tt-component 
of the equation is thus left with three different terms. The metric contributes a term that goes 
like r^, and the gauge and fermion stress tensor contain two different terms: one contributes to 
a r^+^^p^ and another term proportional to qNc2 — J2f goes like ~ r^~^^^p. The only way to 
obtain a non-trivial solution is to switch this latter term off. Therefore we have, 

J ^ n-K {N - l)A'7r 

' = ^t^m^ 2,iv/ ~ 2^' ^^-^^^ 

where I is the regulated size of the time direction. To obtain non-vanishing spinor solution, one 
is forced to take = 0. This relates the index z to the fermion mass, 

z = 2mL - 1, (3.37) 

where we have restored L. With that the Einstein equation is thus reduced simply to an 
algebraic equation constraining pi 

-2N'^h'^ \pi + _ 1)^3 ^ (-3 33^ 

which readily yields 

,2 _ ((mL-l)(2mL- 1)3)1 
l^^l hN-q ' ^^-^^^ 

where we have restored the radius L here. However, pi is the normalization of our plane wave 
solutions. It has to take a fixed value. The normalization of the wave function is determined by 
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d^x^lir)^+ir) = = 1, (3.40) 
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where V is the regulated volume. 

Therefore, the relation (3.39) fixes the value of N that would back-react to give rise to the 
particular geometry and electric fields.i.e. 

{{mL-l){2mL-l)^)-2 _ 1 

hNq V ^ ' 

which gives 

^ = pn = {{mL - l)(2mL - 1)^)^ ~ (3.42) 

where p is the charge density. Einstein equation dictates that ~ 1 which is expected from 
our intuitive argument presented at the beginning of the section. 

Being a constant solution, the fr-component of the fermion energy momentum tensor 

Tl='-^{^pli^+~4,li;'+), (3.43) 

which is a lone term in the Einstein equation also vanishes. Similarly replacing m — >■ —m we 
will get z = —2mL — 1. 

There are a few interesting limits to take. One could for example consider taking the limit 
mL^l such that z ^ Q. Staring at (3.38,3.39) it means a non-trivial fermion solution is only 
recovered for g — t- at the same rate. In this limit however the gauge component At becomes 
a constant, but in fact an infinite constant since C2 is proportional to 1/q. Remarkably this 
is precisely what happens in the gauge-gravitational Chern-Simons analysis when we take the 
limit z ^ and /x^ — > in equation (2.5). Another possibility to consider is that for g = 
to start with, then the homogenous sohition in the Maxwell's equations, namely the log term, 
could have been allowed without entering into Einstein equation through the fermion stress 
tensor. Indeed we checked that this is a solution, provided that m = ^ and 2; = 0, 1. When 
z = 1 the fermion radial wave-function is in fact not a constant but goes like 

^+{r) - r-^ (3.44) 

When one compares with the gauge-gravitational Chern-Simons setting, where one analogously 
takes fip = and allows for the homogenous solution Af ^ C2 + Ci logr, it turns out z = 0,1 
are again solutions, provided that one has to take 

fiG{z - l)z ~ cf (3.45) 

in either case. i.e. /xg curiously approaches infinity. There appears to be a one-to-one cor- 
respondence between the gauge-gravitational Chern-Simons theory and the Einstein-Maxwell- 
Dirac system, as would be expected if Chern-Simons terms should be recovered by integrating 
out fermions. Although our results suggest that the relationship between the value of the 
Chern-Simons couplings and the fermion mass and charge depart from that in flat space. 

3.2 3+1 Dimensions and beyond 

The procedure above can be simply generalized to higher dimensions. Consider for concreteness 
the case D = 3 -I- 1. In this case the metric ansatz we use is the 4-d Schrodinger metric, 

ds^ = L^i-r'^dt^ - Aerdtdx + rdy^ + ^), (3.46) 
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where we will again set L = 1 from now on. Using the same choice of vierbeins as in 2 + 1 with 
the addition of 

= ^/fdy, (3.47) 

where the hat denotes tangent coordinates. The spin-connections will be given again by (3.12), 
with only the new addition of 

a;<^3 ^ -y^Pdy. (3.48) 

Choosing again the Dirac spinor to be 

i,X 

where x = {Ij 2, 3, 4} are spinor polarization indices, 



y^ivMM ~ \ , , (r) r ^±,<^MM - ,/,x,2 J ' K'^-^'J) 



and we take as before ki = 0, i = 1, 2. For clear notations, we will omit the polarization index x 
in the following, only to look explicitly for the non-trivial spinor polarization that would solve 
the system of equations. Assuming only r dependence and going through similar analysis as in 
the previous section leads to very simple results. Starting again with the Maxwell equations to 
obtain constraints on the spinors, the equation for Ar gives 

where as in the previous section wc arc again denoting the components of the tp field operator 
with a hat, to be distinguished from the modes. The simplest solution is again taking, for each 
mode 

This choice also easily ensures the vanishing of the source currents for and Ay. Subsequently 
the Einstein equation and the Maxwell equations would again imply 

.A'p — — — ^* (3.53) 

Fermions taking constant value is forced to be a solution 

for some normalization constants and we arrive at 

A, = c\ log. + 4 + Wl '"^-t(i7|^""'* lQ>. (3.55) 

where 

N+ 7V_ 

1^) =n«un«-,-,io)' {Q\^'Mi\Q) = N±\pi\\ (3.56) 

i i 

and a±^^^^ denotes the creation operators of the fermion modes "0^,0), respectively. N± are the 
fermion numbers for the it modes respectively. As in d = 2 + 1, Einstein equation dictates that 
c\ = and 

{N+IpW + N-\pU')4 = I {j:^Up\\' + I:^-j\P-\'^ ■ (3-57) 
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The two remaining Dirac equations for constant fermions cannot be satisfied simultaneously for 
arbitrary massive fermions when both tp}^ and ■^l are non-vanishing. Explicitly, the remaining 
Dirac equations are given by 

pl{2 + 2m + z) = 0, p\{2 - 2m + z) = 0. (3.58) 

We are thus left with three different possibilities. 

3.2.1 Massless fermions 

For massless fermions, 

z = -2, (3.59) 

and the remaining Einstein equation is 

3(12-(Q|(^it^i + V;yV.|)V)|Q) 



2r2 



0, (3.60) 



which relates iV_ to N+. 



3.2.2 Massive fermions 

For massive fermions, there are two independent sets of solutions. 

,1=0, |-.|. ,2./(9 + 18,„ + 8n.')(^ + l)(2„ + l-F 

g V 4m + 3 

and 

z = 2mL - 2. (3.62) 

This solution applies when m 7^ — |. When m actually takes that value, V'l becomes uncon- 
strained because the Einstein equation vanishes without further constraint on ^1. Similarly 
one could have 



„i _ n N 1.1 |2 - , V(9 - 18m + 8m^)(m - l)(2m-^ 



(3.63) 



where 

z = -2mL - 2 (3.64) 

and m 7^ |. Otherwise A^_|- is unconstrained. 

Here we have an interesting pattern to note. When D = 2+1, the resultant index z = 2mL— 1 
and at D = 3+1 that gives, depending on chirality, ±2mL— 2. It turns out that exact solutions of 
charged constant spinor solutions and Schrodinger metrics can be found in arbitrary dimensions. 
In D = 4 + 1 for example we have 

z = ±2mL - 3. (3.65) 
There is a clear pattern in the constraint on z as we go to higher dimension, namely 

z = ±2mL-{D-2). (3.66) 

The reason it gets shifted in this particular manner is that these constraints follow from requiring 

the coefficient multiplying ^ in the Dirac equation to vanish, for the existence of constant 
spinor solutions. The mass term enters in the same way for arbitrary dimensions, but for each 
additional dimension we get an extra spin-connection component, as illustrated in (3.48), which 
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enters into the Dirac equation, after multiplying by the curved space gamma matrices, simply 
as a constant. Therefore as we move onto higher dimensions the constraint on z get recurrently 
shifted by —1, starting from the value when d = 2+1 ''. The same chiral modes that are obtained 
in the previous section appear also here. However, in dimensions higher than d = 2 + 1, spinor 
dimensions increase such that the constraints from Maxwell's equations become a much weaker 
restriction and it should then be possible generally to switch on momentum in other directions. 
As a result our solution would not be able to capture the back-reaction of a general state where 
modes with non-trivial momenta k are occupied. 



3.3 Can there be an exact Lifshitz solution? 

Given our analysis with the Chern-Simons theory and the form of the probe fermion stress 
tensor, it would be interesting to look for exact Lifshitz solutions. However, in the A;^ = case 
in 2+1 dimensions we can show that this is not possible. Note that in light of our discussion 
of the semi-classical analysis, this will only make sense if the charge q is large compared to the 
mass. The following should only be treated schematically and we leave a more rigorous analysis 
for future work. 

We will again start with the Einstein-Maxwcll-Dirac system in 2+1 dimensions and look for 
an exact solution to the equations of motion in the background given in eq.(2.6). The set of 
equations are 

V^B^^ - mV' = 0, V^F^. = i^, (3.67) 

Ro.^ - IgapR - ^^9a0 - ^(^a7 V " \^^9afl) = ^T^^. (3-68) 

Since we are looking for field configurations depending only on r, the left hand side of the 
equation of motion vanishes trivially. This imposes a chirality constraint on the components of 
the Dirac fermion, namely, 

V'+ = V- ■ (3-69) 

Inserting this condition back into the Dirac equation, we get two equations for ^_|_, which are 
inconsistent unless 

A., = -r'-^At . (3.70) 
Using this condition the Dirac equation can be solved, 

. — z — 2m 

ip+{r) = Vrpir 2 . (3.71) 

We can then substitute this solution and the relation between A^ and At, in the Einstein 
equation. The x-r component of the Einstein equation then implies that ip+ is a real function 
and as a result the integration constant pi is real. Having done this, we can solve the At 
component of the Maxwell equation, which yields 

Atir) = C2 + — — r . 3.72) 

1 — z {l — 2m){2m — z) 

Finally we can solve for the remaining components of the Einstein equations. Although we have 
four equations, only two of them are independent. After having already determined fermion and 
gauge fields, these two equations are purely algebraic in nature and they constrain the remaining 
undetermined integration constants, Ci, C2 and pi and the parameter z. Unfortunately, there is 



'^Note also that in all these solutions T\ — T^^ = 0. 
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no consistent choice of these constants that can satisfy both equations. The same analysis can 
be apphed generally for ^ 0, which suggests that Lifshitz spacctimcs arc not exact solution 
to the equations of motion. This however does not contradict the fact that in the WKB limit 
it becomes an approximate solution. 



4 Discussion 

Motivated by results of [4], we studied a system of self-gravitating charged Dirac fermions and 
showed that non-relativistic Schrodinger metrics exist as exact solutions ^. In 2+1 dimensions, 
this was anticipated by studying the three dimensional Einstcin-Maxwcll-Chcrn-Simons system. 
Although the gauge Chern-Simons term does not couple to the metric, it affects the background 
by modifying the Maxwell equations. We found that in the presence of gauge and/or gravita- 
tional Chern-Simons term, one generically finds Schrodinger spacetime as a suitable background. 

One may be worried that by treating the problem involving fermions semi-classically, we 
could be missing important higher order effects. For starters, the effects of quantum gravity are 
under control as we are implicitly assuming that the gravitational coupling is small. Next, one 
could be worried that since there are a large number of fermions, there could be an enhancement 
of higher loop effects by their presence. However, as demonstrated in a general toy-model in 
Appendix D, as long as the number of modes N excited is large, the tree level contribution 
from contraction with the external state would be boosted by factors of relative to quantum 
loops at each order of perturbation in the boson-fermion interaction. This gives support to our 
semi-classical treatment. 

Since in the flat space, (charged) fermions induce gauge and gravitational Chern-Simons 
terms through parity anomaly^, the result obtained in the Einstein-Maxwell-Chern-Simons was 
a precursor to our expectation of back-reaction of fermions on AdS space. Armed with this result 
we studied fermions in the Schrodinger background in Einstein-Maxwell theory without Chern- 
Simons terms, anticipating that the Schrodinger spacetime is the fully back-reacted geometry 
of the fermions in AdS space. We indeed found a consistent solution to the gravity-MaxwcU- 
fermion system in the Schrodinger spacetime in arbitrary dimension for the critical exponent 
z = ±2mL — (D — 2), where D is the number of bulk spacetime dimensions. 

The emergence of Schrodinger spacetime can also be seen perturbatively. Let us start with 
AdS3 spacetime with metric 

ds^ = L^{-Ardtdx + ^), (4.1) 

and consider gauge field components taking at most constant values. We can then solve the 
Dirac equation for a fermion with charge q and mass m in this background. The solution is 

tp[r,t,x) = r 2 4 exp(— iwi) I ^ I. (4.2) 

We then look at the back-reaction of this solution on the background metric and gauge fields. 
Since we are treating this fermion solution perturbatively, the back-reaction will be tiny (pro- 
portional to h) but indicative of the effect of macroscopic droplet of these fermion on the back- 
ground. With this in mind we evaluate fermion current and energy momentum tensor and feed 



^It is important to note that the fermion gas we are considering is a smooth density distribution, which alleviates 
Coulomb repulsion particularly along the a;-direction, that otherwise would have made the gas unstable. 

^Fermion mass term, in fact, breaks parity symmetry explicitly in 2+1 dimensions. On the other hand in 3+1 
dimensions it does not break parity. 
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it back into equations of motion of the background. Solution to equations with back-reaction of 
fermions generate non-trivial gauge field component At and metric component gu, which take 
the form 



Mr) 




(4.3) 



9tt 



1 / -4g2r-2™-i 



+ ei Inr -|- e2r 



,m-l/2 j^^ 



(4.4) 



(2m -1)2 l^(2m2-3m + l) 



where q are constants of integration and the Cj's are fixed in terms of Cj via Einstein equation. 
In particular, setting ci = sets both = 0. In this limit we are actually reduced to the 
original solution we found previously in (3.10), (3.32) and (3.37). The fermions only appear 
differently because we have picked a different set of vierbeins to suit our purpose here starting 
off from a off-diagonal metric. 

It is interesting to note that for z > the Schrodinger backgrounds are horizon-free and 
as a result do not have macroscopic entropy associated with them. The phenomenon of disap- 
pearance of the horizon due to fcrmion back-reaction is similar to what was proposed in [4]. 
It would be interesting to find a dual description of our solutions. In the 2-|-l dimensions we 
studied the back-reacted geometry supported either by the Chern-Simons terms or by fermions. 
As mentioned earlier, around fiat backgrounds it is known that Chern-Simons terms are induced 
by one loop effects in the fermionic theory [8, 7]. In the context of curved space we have two 
Chern-Simons terms, gauge and gravitational. The Schrodinger metric obtained in the Chern- 
Simons theory has (up to sign) two new exponents apart from the standard one corresponding to 
AdS space; z = and z = {jic T 1)/(2mg) (see (2.5)). On the other hand, in the fermionic 

theory wc have only one exponent (up to a sign of the fermion mass term). It is tempting to 
conjecture that even in curved spacctimes, the fermions would induce Chern-Simons terms. It 
would be interesting to directly relate these two theories and compare their resulting exponent 
z. 

Let us make a brief comment on the Maxwell-Chern-Simons theory. Around flat spacctime 
it is possible to arrange the coefficient of gauge Chern-Simons term and the number of fcrmion 
species to switch off the gravitational Chern-Simons term[9]. It turns out that this system has 
already been studied in [21, 22] and [23]. In [23], it was proved that there exist no black hole 
solutions in this theory. The existence of Chern-Simons term in 2+1 dimensions implies the 
presence of a chiral anomaly in the 1-1-1 dimensional theory [24]. The induced chiral anomaly 
in the boundary theory is canceled by the chiral current in the boundary theory. In [22], it 
was found (after transforming to our conventions), that regular black hole solutions only exist 
if hfI^G > 1/6. This is reminiscent of the condition (C.ll). 

It has been shown in [4] that Lifshitz metrics are solutions when a WKB approximation is 
made. We leave it as an open problem to examine if they are also exact solutions. In light of 
our analysis using the Chern-Simons theory and the self-gravitating fermions, this is unlikely 
to happen in 2-1-1 dimensions. However in higher dimensions, there are more ways to solve 
the constraint j.,- = than in 2+1 dimensions. So it could very well be that both Lifshitz and 
Schrodinger metrics exist as exact solutions. It will be exciting to find new backgrounds from 
self-gravitating fermions. Of course in general, this is a very difficult problem. 
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A Conventions 



We use the following conventions for three and four dimensions. The three dimensional action 
with Chern-Simon terms is 



2ipl 



g{R+^- ^Fi^.F'^'') + Lmcs + Lgcs 



where the Maxwell Chern Simons and gravitational Chern Simons are given by 

Lgcs = i^ce'^^^n^idA + 1^^,^^^ ■ 



The four dimensional action is 



2^2 



-g{R + 



L2 



The Einstein equation is 



(A.l) 

(A.2) 

(A.3) 
(A.4) 

(A.5) 

(A.6) 



where A is the cosmological constant in appropriate dimensions and r^,^ is the matter stress 
energy tensor. The matter stress energy tensor with this choice of conventions is defined as 



5S 



The generic asymptotically AdS metric is written as 



r2 /(r) L2 



+ -T^dxj . 



(A.7) 



(A.8) 



In absence of Chern-Simons terms in three dimensions it can be shown that (A.8) solves equa- 
tions of motion to give charged BTZ black hole solution[25, 26, 27] with 



2 r2 



/(r) = l- 



<fL^ q'L 
4^2 



In—, At = gin—, 
2r^ r* To 



(A.9) 



where ro is the location of outer horizon. The condition for existence of two horizons is g^L^ > 
4rl and the extremality condition \s = qL/2. 

In four dimensions, (A.8) solves the equations of motion to give charged Reissner-Nordstrom 
black hole solution, with behaviour of /(r) and At{r) given by 



, , M q^L"^ , , , ,1 1 , 
/ r = 1 - — + Mr) = qi 



(A.IO) 
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In order to have double horizon 27M^ > Aq^L^ and the equaUty corresponds to extremal 
solution. 

The near horizon metric in the extremal limit is given by AdS2 x Af„, where M„ is an n- 
dimensional manifold, which in our case is either i?" or S". Our convention for AdS2 metric in 
Poincare coordinates is 

^^4^52 = -T^dt^ + ^dr^, (A.ll) 



where L is the radius of AdS2. In 2+1 dimensions, according to our conventions, 

At = V2l, L = A, (A.12) 

while in 3+1 dimensions, 

At = V2l, L = A. (A.13) 

B Spin connection and the Dirac Equation 

Consider a general metric of the form 

ds^ = -gtt{r)df + grr{r)dr'^ + gxx{r)dx^ + gii{r)dx'dx' + 2gtx{r)dtdx , (B.l) 

with all metric components having only radial dependence. We will use /i, u, - ■ ■ to denote 
curved spacetime indices and to denote curved spatial indices. We will reserve index 

'd' for denoting flat coordinate related to curved index r by vielbein. Flat space indices for 
boundary coordinates are denoted by a, 6, • • •. 

In this case, following choice of vielbeins is suitable for both diagonal metrics as well as 
Schrodinger metrics 

e° = ^/gttdt - -^^dx, e'^ = Jg^dr 



1 _ \l9tt9xx + 9l ^^^ ^„ ^ ^idx\ a = 2, 3, • • • . (B.2) 



/Qtt 

It is straightforward to derive components of spin connections from them and they are given by 
Wdo - o , dt - - dx, 

^ V 9tt9rr ^ V 9tt9rr 9tt 

^di = ^ 2^{-dctg J^t t^^"^^^^'''' ^ 3tx9'tx{r))dx + {9tt9'txiT) - 9tx9'tt{r))dt] , (B.4) 

uJdi = /'^ dx\ (B.5) 

where primes one the metric components are derivatives with respect to their argument. 

Let us now look at the Dirac equation in different backgrounds. The Dirac equation for a 
particle with mass m and charge e takes the form 

{p - m)V' = 0, (B.6) 

where 

p = V^e^id^ + ^uj^abT'^T' + iqA^). (B.7) 
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To solve the Dirac equation we will use the following F-matrix convention for 2+1 dimensions, 

T^ = a3, r0 = ic72, r'^ = ai. (B.8) 
Our 3+1 dimensional F-matrix convention is 

Where 7 is 2 x 2 identity matrix and cr" are Pauli matrices. 

C Stress tensor of fermionic probes in AdS 

We will be considering fermions in curved space. Wc shall begin presenting-^'' explicitly the 
fermionic contribution to the stress tensor [19]. The quadratic action for a charged massive 
Dirac spinor is 

^Dira. = -ir / d'^+'xV^ii^r'^e^D^i; - m^V) (CI) 

where 



Vi = V^FO, F„;. = --[r„,F,] (C.2) 
The corresponding stress tensor is 

= \ + iV;^(^F,)V') , (C.3) 
and symmetrization is denoted by V^^Vu) ~ ^(^^K^ + 

C.l AdS2 X R 

The background is, following our conventions, 

ds\dS2xR = -T^dt^ + -^dr^ + dx^. At = (C.4) 



12 ' ■ 

and the only non-vanishing spin connection is. 



mo = ^dt. (C.5) 

Taking a Dirac spinor of the form 

and assuming only r dependence, the Dirac equation is given explicitly by 

2rV'± + 2iqV2ip± + (1 + 2mL)V'=F = 0. (C.7) 



^"Note that [19] use the mostly negative convention while we follow the more standard GR convention of mostly 
positive. 
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The solution is simply 

for some constants p^. Switching on one of the two solutions leads to 



Ti 



f _ \pf\'^V2m{mL^5) ^^^s 



tt 



Note the appearance of T^j components. The rr component contributes only from cross terms 
when both solutions are switched on. Note that the Ricci tensor for the background metric is 
simply 

r2 1 
Rtt = Rrr = (C.IO) 

Therefore, in the near horizon limit r — > 0, the solution scaling with r~^2+'') would eventually 
become dominant over the background such that the probe approximation breaks down. For 
the other solution that scales as r^^a"*^) however, the probe approximation is good all the way 
to the horizon if (5 > 1, which gives 

rr? > 2q^ + 1. (C.ll) 

When the mass is smaller than 2q^ + 1, back reaction in the near horizon region cannot be 
ignored, with the likely result of destroying the horizon altogether. Therefore in the large mass 
and charge limit we recover the assertion in [4] that roughly horizon is destroyed whenever eq > 
m, where e = \/2 is the background electric field. Note also that these solutions independently 
satisfy 

Tl - T; = _\pt\'V2rnimLTS)^_,^,s < ^ (^.12) 
Lq 

C.2 AdS2 X 

Let us repeat the exercise in the previous section in AdS2 x E?. A static solution of a 4- 
component Dirac fermion would be given by 

^ = V;+r-i+'^ + V-r-"^"^ 5 = L^Jm? - 2q^, (C.13) 

where 



*l = ±fc(..*i), (C.14) 

where -^j^ are two component constant spinors and CTj are the Pauli matrices. The non-vanishing 
contribution of ^+ to the energy-momentum tensor is given by 

T,^ = -^^a^Vi (C.17) 

(C.18) 



m] 
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As in the previous section, the energy-momentum clearly illustrates again that back-reaction is 
severe in the near horizon region generally for \/29 > '"^ Note that here 

C.3 AdSj X 52 

An entirely parallel story goes through if one were to replace in the previous sub-section by 

S^. Namely 

^ = ^^r-5+'5 + ^_r-5-'5, § = L.Jm'^ - 2q'^, (C.20) 

where, however even the lowest mode contains non-trivial angular dependence. 
The energy-momentum tensor here is similarly. 



Tf ^ csc.,lt,.:!f!!(^!(l^^ 
Tf = csc^^Vs^r'"^^'^"' - ^ ^ '^'^^ (^•^^) 

^"^ 8qr " 2V2L^ ^ ' ' 

(C.25) 

The null energy condition is the same as eq. (C.19). 



D A toy model justification for our semi-classical 
treatment 

Our discussion is largely based on the assumption that in the presence of a large fermi-gas, one 
can treat the background classically, while the fermions propagate in this mean-field background 
quantum-mechanically. We consider here a toy model where fermions are coupled to bosons. 
In the case where the fermion fields appear only up to quadratic order, we can obtain the 1- 
loop effective action of the bosons after integrating out the fermions in an A?^-fermion state. 
It turns out, in the large iV-limit and to leading order in the boson- fermion coupling g, the 
effective equations of motion of the bosons are equivalent to the original one, except the fermion 
operators are replaced by their expectation values evaluated on the A'^-fermion state. The toy 
model gives us more confidence in the methods we advocate in this paper. It is also noteworthy 
that the approach of treating back-reacting fermions on a mean-field background and obtaining 
self-consistent set of solutions to all the equations of motion with an explicit built-in Fermi- 
surface, is a well-known strategy in nuclear physics[15], which is equivalent in various limits to 
the Hartree-Fock approximation. 

Our discussion here will follow closely [28]. The (Euclidean) action of our toy-model at 
temperature T = ^~^in (0-1-1) dimension is as follows: 

1 ff^ 

S = Sb + Sf + Sint, Sb = - dT4>^ + ml4>'^, 

^ Jo 
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Sf + Sint = dTj2 blh + {nii - Hi)b\hi - g<ph\bi, (D.l) 



where we have introduced many species of fermionic fields hi and a corresponding chemical 
potential for each species, which are to be determined from some given occupation of the 
fermionic states. The partition function is given by 



Z = J D[cP]l[D[zbl]D[h]eM^)- (D-2) 



Since the action is quadratic in the fermions, the path-integral can be readily done. Expressing 
the fields as 

bi = Y. = E ^"'"""'^n (D.3) 
n n 

The cOn here are thus Matsubara frequencies given by 

Loi = (2n + l)7rr, to';, = 2nnT, (D.4) 
and Sf + Sint can be re- written as 

Sf + Sint = ^^im^mnhn, T^mn = -i/3i{-iuin + l^i - mi)^mn - 94>m-n)- (D-5) 

i,m,n 

Integrating out the fermions one is then left with 
Z = j Z^[0]exp(-|)det(|) = det(^Ha;„+/x,-m,)) j D[ct>] exp(-|) det(l- _ ^ ). 

To lowest order in the coupling the fermion contribution to the free energy evaluates to [28], 
using various summation formula, 

p. 

= - In = - 2J ln[^(-za;„ + iJ.i- rrii)] 

n.i 

^ i n 



where we have dropped the contribution of vacuum energy which has to be regulated. We can 
determine the chemical potentials by 

which is precisely the Fermi-distribution. The ground state 

iQ)=n^^io)' (D-9) 

i 

corresponding to a system of N fermions with fermi-energy Ep can thus be obtained by picking 

^^^ = Ep. (D.IO) 



22 



In which case the Fermi-distribution behaves Uke a step function, giving occupation ATj = 1 for 
all single particle states below Ef. Now returning to the full path integral, we have 

Z =1 D[</>]exp(-^)exp(lndet(|)) 

~ J D[^] exp(-y) exp(^^mi + ln(l + e-^('"*-^/-s^o)) - In^), (D.ll) 

i 

where in the second line we have taken the leading g approximation, so that in the determinant 
of a matrix of the form 

(1 +#o) 94>i ■■■ \ 
gcp-i (l + #o) ••• 1 (D.12) 

the leading g dependence came from the trace, which include only the diagonal components 
with 4>(m-m,=o) ■ The off-diagonal components can only arise in higher g corrections. 

Since in the zero-temperature limit /3 — >■ oo the factor e~^^'^^''^f^ approaches infinity for 
rrii < Ef and zero otherwise, we have 



Z ^yz)[0]exp(-^)exp( ^ _f3^^. _ Ef - gcj^o)) 

i=filled states 

~ I D[^] exp(-^) eM-P{NgM + •.•), (D-13) 

where the vacuum energy that is supposed to be regulated, and other pieces independent of (p are 

denoted in Note that to leading order in g the effective action of the bosons after integrating 

out the fermions, is precisely as if we replace the fermion operators by their expectation values 
in the A/^-fermion state. For back-reaction on the bosons to be important however, we need 

Ng ~ i. (D.14) 

Note also that the approximation we are taking in the partition function is equivalent to, 
from the perspective of perturbation theory via Feynman diagrams, only the tree diagrams i.e. 
The perturbation operators are only contracting with external states. Consider for example 
the leading order g correction to any expectation values of some arbitrary operators O in the 
presence of the interaction terms, 

{Q\0 1 4>Y.ih\Q), (D.15) 

■' i 

we see that contraction of the bi's with the external creation operators in \Q) result in a factor 
of N, whereas contraction among the 6j's, corresponding to loop diagrams, is only of order 1 
relatively. Similarly at order g"^, we have 

^// T.<l>mOiblh)ib]bj)\Q), (D.16) 

where contraction of the 6's with the external state would again give rise to an extra factor of 
roughly N{N — 1) ~ A^^ relative to contraction among the b and 6^.^^ In the limit where N is 
large and g is small, the loop diagrams would be suppressed relative to tree diagrams. 



""^^Note that the combinatoric factor of 2! would appear in either internal contractions or contraction with the external 
state. We have also made used of the fact that states with different excitations are orthogonal to each other. Therefore 
having chosen a particular contraction for the &j's on the contraction of bj with {Q\ is completely determined. 
This explains why there is only a factor of N{N — 1), rather than the square of it. 
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To summarize, in the limit where Ng ~ , the back-reaction of the fer mi-gas on the 
bosons become significant, but at the same time we have the quantum loops under control 
because g is small. Particularly, the leading g effect on the bosons can be captured by replacing 
all the fermionic operators by their expectation values evaluated on the filled state \Q) in the 
equations of motion of the bosons. In fact the coupling to the bosonic field could have 

been generalized to some arbitrary one of the form gV{(t>)b^b without altering the conclusion, 
as long as g is small. The toy model appears to justify the approach taken in this paper in 
including fermionic back-reactions on a bosonic background. 
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